Non-minimal Split Supersymmetry 



S. V. Demidov ' 1 , D. S. Gorbunov a > 2 
a Institute for Nuclear Research of the Russian Academy of Sciences, 
60th October Anniversary prospect la, Moscow 117312, Russia 



Abstract 

We present an extension of the minimal split supersymmetry model, which is 
capable of explaining the baryon asymmetry of the Universe. Instead of MSSM 
we start from NMSSM and split its spectrum in such a way that the low energy 
theory contains neutral particles, in addition to the content of minimal split su- 
persymmetry. They trigger the strongly first order electroweak phase transition 
(EWPT) and provide an additional source of CP-violation. In this model, we esti- 
mate the amount of the baryon asymmetry produced during EWPT, using WKB 
approximation for CP-violating sources in diffusion equations. We also examine the 
contribution of CP-violating interactions to the electron and neutron electric dipole 
moments and estimate the production of the neutralino dark matter. We find that 
both phenomenological and cosmological requirements can be fulfilled in this model. 
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1 Introduction 



In spite of approximate symmetry between particles and antiparticles, our visible Universe 
is asymmetric in baryons. An explanation of this fact ought to be addressed in any viable 
theory. Quantitatively, the measurements of the anisotropy of the cosmic microwave 
background constrain the baryon-to-photon ratio [1] to be in the range 

6.1 x 1(T 10 < — < 6.9 x KT 10 . (1) 

Three necessary conditions (Sakharov conditions) must be fulfilled in the early Universe 
to produce the baryon asymmetry [2] : baryon number violation, C- and CP- violation and 
departure from thermal equilibrium. 

Electroweak baryogenesis [3] is one of the most interesting scenarios of baryon asym- 
metry generation. If the electroweak phase transition (EWPT) is of the first order, it 
provides departure from thermal equilibrium in cosmic plasma, as it proceeds through 
nucleation of bubbles of the broken phase, their subsequent growth and percolation. The 
baryon asymmetry generated during the phase transition can be washed out by rapid 
sphaleron processes in the broken phase [4]. In the Standard Model, the condition 

V 7f > 1-1 (2) 

1 c 

(where T c is the critical temperature and v c is the Higgs expectation value at this tempera- 
ture) would ensure that the sphaleron transitions are suppressed inside the bubbles of the 
new phase [5] and the baryon asymmetry survives. Nonperturbative analysis [6] revealed 
that, with the present experimental bound on the Higgs boson mass, the electroweak 
phase transition in the Standard Model not only does not yield (2), but is definitely not 
of the first order. Moreover, in spite of the presence of CP-violating phase in CKM ma- 
trix, its contribution to the production of the baryon asymmetry is too small to explain 
the observed value (1). The Standard Model fails to support electroweak baryogenesis. 

Supersymmetry is one of the most attractive ways beyond the framework of the Stan- 
dard Model. One recalls the cancellation of quadratic divergences and gauge coupling 
unification among the main reasons for interest in supersymmetric theories. Minimal 
Supersymmetric Standard Model (MSSM) and its extensions suggest the lightest super- 
symmetric particle (LSP) as a natural dark matter candidate and provide mechanisms to 
generate the baryon asymmetry (for recent studies see Refs. [7], [8], [9], [10] and references 
therein). 

Recently, motivated by landscape of vacua in string theory and cosmological constant 
problem, models with split supersymmetry have been proposed [11,12]. The spectrum 
of these models is governed by two scales. While the electroweak scale m ew determines 
the masses of the Standard Model particles and additional fermionic degrees of freedom 
(higgsinos and gauginos) the other MSSM particles (scalars) have masses of order of a 
splitting scale m S) which generally may be in a wide range of 10 4 — 10 15 GeV. Such a 
pattern ensures the gauge coupling unification, so that split SUSY can be incorporated 
into a GUT. 
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Split SUSY enables one to avoid phenomenological difficulties with CP- and flavor 
violation and also provides natural dark matter candidates. But from the point of view of 
the electroweak baryogenesis 3 , split super symmetry, in its minimal version, shares some 
disadvantages of SM: the electroweak phase transition remains too weak, because the 
additional, weakly coupled fermions do not strengthen EWPT. It is of interest to find a 
viable scenario, which would be compatible with general split supersymmetry framework 
and would provide conditions required for effective electroweak baryogenesis. One way, 
proposed in Ref. [14], is to spoil the supersymmetric relations between the gauge and 
Yukawa gaugino-higgsino couplings. For rather large Yukawas, the charginos and neu- 
tralinos are light in the symmetric phase but heavy in the broken phase, so they decouple 
from the cosmic plasma at the moment of the phase transition. During the EWPT the 
particle species to be decoupled transfer their entropy to the cosmic plasma, that results 
in the increase of the temperature and the delay of the transition, providing a possibility 
to satisfy the inequality (2). 

In this paper we pursue another approach and propose a model which keeps all phe- 
nomenologically interesting properties of the minimal split supersymmetry and at the 
same time is capable of producing the baryon asymmetry at EWPT. Instead of using 
MSSM as a starting point, we begin with Non-Minimal Supersymmetric Standard Model 
(NMSSM) and split its spectrum in such a way as to provide the low energy theory with 
all features required for successful electroweak baryogenesis. Indeed, additional neutral 
scalar fields, which are left at low energies after the spectrum splitting, change zero tem- 
perature potential and can make EWPT stronger. The model contains new sources of 
CP-violation as well. Apart from their crucial role in the generation of the required 
amount of the baryon asymmetry (1), they also contribute to the electric dipole moments 
(EDMs) of electron and neutron at the level detectable in the near future experiments. 
Also, we study the possibility that the lightest neutralino is dark matter particle. We 
find that along with gaugino- and higgsino-like candidates there is a region of parameter 
space where dark matter is mostly singlino. 

The outline of this paper is as follows. Sec. 2 contains the description of our model 
and its main features. In Sec. 3 we study the strength of the electroweak phase tran- 
sition and show that EWPT can be the strongly first order. In Sec. 4 we describe the 
baryon asymmetry calculation procedure and present numerical results. In Sec. 5 we es- 
timate the values of electron and neutron EDMs in our model. Neutralino dark matter 
phenomenology is discussed in Sec. 6 and Sec. 7 contains our conclusions. 

2 The model 

We begin with the most general NMSSM. The relevant for our study part of the super- 
potential is [15] 

W = \NH u eH d + -kN 3 + fiH u eH d + rN. (3) 
o 

3 For a discussion of prospects of the Affleck-Dine mechanism of baryogenesis in split SUSY see Ref. [13]. 
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Here 4 H u and H d are the Higgs doublets, N is a chiral superfield, which is a singlet with 
respect to the SM gauge group and e is antisymmetric 2x2 matrix, e 12 = 1. Soft SUSY 
breaking terms for this model read 

V soft = (\A x NH u eH d + l -kA k N* + fiBH u eH d + A r N + h.c)j (4) 

+m 2 u HtH u + m 2 d H d H d + m 2 N \N\ 2 . 

Electroweak baryogenesis in NMSSM has been explored in Refs. [8] in the context of 
low energy supersymmetry breaking. We consider this model in the framework of split 
SUSY. For this purpose some parameters need to be fine-tuned, like in the minimal split 
SUSY, in such a way that the particle spectrum is split into two parts. To strengthen 
the electroweak phase transition, some scalars can be recruited [16]. To preserve gauge 
coupling unification, these particles should not contribute to the beta functions (or their 
contributions should be canceled) at least at the leading order. So in the minimal case, 
the low energy theory can contain (besides the usual split SUSY particles) singlets with 
respect to the SM gauge group. 

We consider the theory in which the only new source of explicit CP-violation is //- 
parameter which for concreteness we take to be exactly imaginary. All other parameters 
(including 5/x-term) of the theory (3), (4) are taken to be real (except in the analysis of 
the electron and neutron EDMs). This is a purely technical assumption which we make 
in order that the parameter space to be tractable. 

We examine splitting in the {H u , H d) N) sector only, because decoupling of squarks 
and sleptons is supposed to be similar to the case of the minimal split SUSY. The tree 
level scalar potential of NMSSM is 



V = V D + V F + Vsoft , (5) 



where 



V D = 9 - (H\o a H d + Hla a H u ) + 9 — {\H d \ 2 - \H U \ 2 )\ (6) 
V F = \\H u eH d + kN 2 + r\ 2 + \XN + fi\ 2 (hIh u + HlH d ^j , (7) 

and the Higgs doublets are H u = (if+,iJ°) T and H d = (H d ,H d ) T . Let us choose the 
vacuum expectation values for neutral scalars as follows, 



(H° d ) = -j=v d , (H° u ) = ^=w u e*-, (N) = -j. (v s + iv P ) = ^4 + v 2 P e^ . 

Here we have taken into account that the effective potential depends only on the sum 
of the phases of vev's of the Higgs doublets and one of these phases can be set equal to 
zero by a gauge transformation. The stationarity condition with respect to the phase 4>h, 
dV/d(pH = 0, gives the following relation 

f (v 2 s + v 2 P ) sin {4> H - 20 s ) + y/v 2 s + v% sin (0 5 + 4>h) 
Sin ^ = WTXr • (8) 



4 We denote superfields and their scalar components by letters with hat and without hat, respectively. 
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To split the spectrum of the theory, we take soft SUSY breaking parameters Bfi, m 2 and 
m 2 , of the order of the splitting scale m 2 and all other parameters in the scalar sector 
(including all vev's) of the order of the electroweak scale m ew . Equation (8) implies 
the estimate \4>h\ ~ (memM) 2 and therefore this phase can be safely neglected in the 
following analysis. 

We find that vev's of the scalar fields at the minimum of the potential (5) must satisfy 
the following conditions 

g 2 X 2 X 2 ( A \ 2 

m 2 u = —v 2 cos2p-—v 2 cos 2 p-YV 2 s -{\^\ + ^vp\ + m 2 A cos 2 (3, (9) 

g 2 X 2 X 2 ( X \ 2 

m 2 d = -—v 2 cos2P-—v 2 sm 2 p~—v 2 s -{\fi\ + —v P ) + m 2 A sin 2 /3,(10) 

1 A 2 XA v 2 

m 2 N + 2kr = -kXv 2 sin 2(3 - k 2 (vl + vp) v 2 + — sin 2(5 (11) 

2 li ^ 2 2^/2 v s 

kAfc ^ v f>\ 



v s - 



m 



N 



V2 V v s J v s ' 
2kr = -hxv 2 sm2(3-k 2 (v 2 s + v 2 P )-^v 2 -^^ + V2kA k v s v P , (12) 



where we use the notations tan/5 = v u /vd, v 2 = v 2 + v\, g 2 = g 2 + g' 2 and 
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(^Xr + ^-(v 2 s - v 2 P ) + -^=XA x vs + pB^ . 



sin 2(3 

We expand the scalar fields about their vev's as follows, 

Rei^ = -L(v d + h° cos (3 - H° sin (3) , \mH° d = (A sin (3 + G° cos (3) , 
v 2 v 2 

Re# ° = (u« + h° sin /3 + H° cos /3) , Im#° = (A cos (3 - G° sin /3) , 
v 2 ' v 2 

# d = G^cos^ + iZ-sin^, = # + cos/3 - G + sin/3, 

N= (v s + S + i(v P + P)) . 
v2 

Here and G° are Goldstone modes which are eaten, due to the Higgs mechanism, by 
W ± - and Z°-bosons, respectively, h° and H° are the neutral scalar Higgs bosons, A is 
the Higgs pseudoscalar, S and P are singlet scalar and pseudoscalar fields. 

Let us choose the basis in the space of neutral scalars as (h°, H°, A , S, P). By making 
use of the conditions (9) - (12), the elements of 5 x 5 squared mass matrix for neutral 
scalars read as follows, 

M 2 n = £^ cos 2 2(3+— sin 2 2(3, (13) 

M 2 2 = m 2 A -*f sin 2 2(3 sin 2 2(3, (14) 
M 3 2 3 = m 2 A , (15) 

M 2 4 = 2k 2 v 2 s + ^-sm2f3 + ^(v s + V -i)-V2^ (16) 
2V2v s x/2 V v sJ v s 



V2 v P 

M 2 2 = M 2 21 =(^- 9 -^jv 2 sm4(3, (18) 

M 2 3 = M^ = -kXv s v P + ^v P , (19) 

M 2 4 = M\ x = -k\v s v sin 2(3 + X 2 v s v - sin 2(3, (20) 

v2 

M 2 5 = M 5 2 4 = kXvpv sin 2/3 + v^Av + -^p) » ( 21 ) 

M 2 2 3 = M 3 2 2 = 0, (22) 

M 2 4 = M 2 2 = -fcAv s vcos2/?- ^—^ cos 2/3, (23) 

v2 

M 2 5 = M 2 2 = A;Aw P w cos 2/3, (24) 

M 3 2 4 = M 2 3 = -kXvvp, (25) 

M 3 2 5 = M 2 3 = -kXvv s + ^v, (26) 

M 2 5 = M 2 4 = 2k 2 v s v P - V2kA k v P . (27) 

For the squared masses of the charged Higgs bosons one has rn 2 H± = m 2 A — pp- + 
From the expressions for the scalar mass matrices it can be easily seen that, like in the 
minimal split SUSY, one can split the particle spectrum by taking \xB parameter of order 
of the squared splitting scale m 2 . This means that the charged higgses, one scalar H° 
and one pseudoscalar A from the Higgs sector become heavy. It is straightforward to 
check that the low energy effective Lagrangian is obtained by the same substitution for 
the Higgs doublets as in the minimal split SUSY [12] 

H U ^H sin (3, H d -> eH* cos (3, (28) 

where H is the Standard Model Higgs doublet. After the replacement (28) and including 
interactions with fermions, we arrive at the Lagrangian for light degrees of freedom 

C = Cy + jCy- 

It consists of two parts: the scalar potential 

—2 \ 

-C v = 9 — cos 2 2(3 (H^H) 2 + \r + kN 2 -- sin2(3H*H\ 2 + \XN + n\ 2 H^H 

8 ' 2 



(29) 



+ (--Ax sin 2(3NH ] H - ^ sin 2(3H ] H + -kA k N 3 + A r N + h. 
V 2 2 3 

+ (m 2 sin 2 (3 + m 2 d cos 2 (3) H ] H + m 2 N \N\ 2 
and the Yukawa interactions 

-C Y = -XNH u eH d - X sin(3H T e (^H d nj + A cos/3 (nH^j H* - kNhh + h.c, (30) 
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where H u and H d are higgsinos and n is singlino, which is the fermionic component of the 
singlet chiral superfield N . In the above expressions we omitted terms with the Yukawa 
couplings of quarks to the Higgs bosons, gaugino interactions and mass terms, which are 
the same as the ones in Ref. [12]. For this Lagrangian to describe the low energy physics, 
all terms are to be of the order of electroweak scale. Inspecting the Higgs mass term 
in (29) (see also eq. (33) below) and eqs. (9), (10), one can find that a strong cancellation 
in the combination m 2 tan (3 + m 2 d cot (3 — 2Bfi has to take place. 

A special feature of our model, as compared to the minimal split SUSY, is that after 
splitting there remain relatively light singlet fields: complex scalar iV and Majorana 
fermion n. 

It is often assumed that NMSSM possesses Z3-symmetry, which implies fi = r = A r = 
B/j, = 0. It was argued in Ref. [17], that the splitting of Z 3 -symmetric theory, which 
leads to the same particle content as the one described above, results in the relations 
A ~ (m ew /m s ) 2 and (N) ~ m 2 /m ew ; therefore, the cubic term in the scalar effective 
potential behaves effectively like quadratic during the electroweak phase transition. This 
means that EWPT cannot be strong enough in that case and electroweak baryogenesis 
does not work. Thus we do not impose the Z 3 -symmetry and take all parameters in (3) 
and (4) to be non-zero. 

The Lagrangian (29), (30) describes interactions at the splitting scale m s , which is 
taken below to be 10 9 GeV, if not stated otherwise. To obtain the low energy theory, the 
couplings in (29), (30) should be changed according to the renormalization group equations 
(RGE). The set of RGE for dimensionless couplings is presented in Appendix A. 

Below the splitting scale m s , the theory is described by the following Lagrangian 

-C y = -m 2 H ] H + - (H ] H) 2 + iAiH^H (N - N*) + A 2 H^H (N + N*) + n x \N\ 2 H ] H 
z\ 

-k 2 H ] H (N 2 + N* 2 ) + m 2 N \N\ 2 + Ajv|iV 2 | 2 + \h (iV 3 + N* 3 ) + A r (N + N*) (31) 

(Tfi^ 1 — 71 

— N 2 + -A 3 N 2 N* + + -N 3 N* + h.c. 
2 2 s 6 

and 

-C Y = -^W a W a + -^-BB + (/i + kN) H^eH d - kNnn 

+if t (^=~g u a a W a + -±=g' u B - \ u h^j H u (32) 

+H T e ^-i^r + ^=~g' d B - \ d n^j H d + h.c, 

where we have also written explicitly the interactions between the Higgs bosons, hig- 
gsinos and gauginos. The terms, analogous to ones with couplings A 3 ,^,r] in the last 
line in eq. (31), are absent in eq. (29), but they are generated below m s by quantum 
corrections coming from one- loop diagrams of the types presented in Fig. 1. Comparing 
the Lagrangians (29), (31) and (30), (32), we read off the following matching conditions 
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Figure 1: Examples of the Feynman diagrams, which generate terms in the last line in 
eq. (31). 



between the coupling constants, which are valid at the splitting scale m s , 

m 2 = - (ml sin 2 (3 + m 2 d cos 2 (3 - /j,B sin 2(3 + \fi\ 2 - Ar sin 2/3) , 
A 1 = A/i, A 2 = —\A X sin (3 cos (3, 
Ki = \ 2 , k 2 = — \k sin /3 cos /9, \n = k 2 , 

■^k ~~ ^■^■k'l ~~ ' ? A r — — ; 7TI j\j — — 777/ ) 

o 2 A 2 
A = ^- cos 2 2/5 + — sin 2 2/?, « = A, 
4 2 

A u = A cos /3, Ad = —A sin /3, 

i 3 = ^ = e = 

and, like in the minimal split SUSY, 

g u = gsin(3, g d = gcos(3, g' u = g' sin/3, g' d = g'cos(3. 



(33) 

(34) 

(35) 

(36) 
(37) 

(38) 



The aforementioned cancellation in (33) results in m 2 ~ m^. The matching conditions 
provide the initial values for RGE. We use 2-loop RGE for gauge couplings and 1-loop RGE 
for other dimensionless couplings, neglecting threshold effects. The coupling constants, 
additional to the minimal split SUSY, begin to contribute to the running of the gauge 
couplings only at 2-loop level and do not spoil their unification. 

In the next section we investigate EWPT in this model and study the properties of 
one- loop effective potential. Its zero temperature part [21] is written in DR scheme as 
follows, 



64vr 2 



n-imf 



log-^--- 
q z 2 



(39) 



where the sum runs over all particle species of field-dependent mass with rij degrees of 
freedom. Upper and lower signs correspond to the boson and fermion cases, respectively. 
The renormalization scale q is chosen to be 100 GeV. We include loop contributions from 
top quark, gauge and Higgs boson, singlet (pseudo)scalars. Let us single out the terms in 
V tree = — £y, which are quadratic in vev's v, v$ and vp, e.g., 



(40) 



8 



where the last term denotes the rest of the tree level potential. Imposing the conditions 

-'" 2 + \l + v ' (1) ) = ■ (41) 
mJ + * 2 + ilk (vil + v " 1,) = m ' s - * 2 + v P ik {v >™ + yW) = ' (42) 

one ensures that minimum remains at (v,vs,vp). 

Let us describe the part of the parameter space, relevant for our study. There are 
three dimensionless parameters tan/3, A, k, which a priori take any values compatible 
with the weak coupling regime in both high-energy and low-energy theories. In addition, 
the low-energy spectrum has to be phenomenologically viable. To reduce the parameter 
space of the model, in what follows we fix these dimensionless parameters at the splitting 
scale as (unless stated otherwise), 

tan/3 = 10, A = 0.6, k = -0.5. (43) 

There is nothing special in this set, except of the choice of rather large value of tan (3. As 
it will be argued in Sec. 4, if tan/3 ~ 1, produced amount of the baryon asymmetry is 
zero within the semiclassical approximation for CP-violating sources. Quantitatively, the 
numbers we obtain (the exact amount of the baryon asymmetry, predictions for EDMs, 
etc.) are stable with respect to small variations of these parameters about the point (43). 

The dimensionful parameters are taken at the electroweak scale. So, we do not assume 
universal boundary conditions for soft supersymmetry breaking terms. Below we vary 
some of the parameters keeping several relations between them to simplify the analysis. 
For m 2 , m 2 N and m 2 we substitute the corresponding expressions in terms of v, Vs and 
vp, which follow from eqs. (41), (42). Since we take \i to be imaginary and CP is violated 
in the Higgs sector, two scalars and one pseudoscalar generally mix. We restrict our 
considerations (the only reason is to reduce the number of free parameters) to the case 
where mixing between the physical Higgs boson and two other (pseudo) scalars is absent; 
in fact this can be achieved by tuning trilinear constants A X) A 2 . Nonzero A 3 is generated 
by radiative corrections, hence generically it is small and is not very important, so we set it 
equal to zero. We also choose A r = for concreteness. The free dimensionful parameters 
we are left with are vev's of singlet scalars v s , v P) gaugino masses M 1: M 2 , M 3 and 
trilinear coupling constant A^. 

We use the RGE as follows. We take known values of the gauge and top Yukawa 
couplings derived from the following relations 5 : a s (Mz) = 0.117, yt{m t ) = 0.95, Mz = 
91.19 GeV, Mw = 80.42 GeV. These four couplings run up to the scale /i = m s according 
to truncated RGE, in which we leave only y t , % — 1,2, 3. At the splitting scale we take 
the set of other parameters in (31) and (32), obeying initial conditions (43) and matching 
conditions (35-38). Then all couplings run down to the electroweak scale in accordance 
with the full set of RGE. Because of additional coupling constants in RGE, the procedure 
described above does not give the correct low energy pattern of gauge and top Yukawa 

5 Thc adopted value of y t corresponds to the top quark mass 171.6 GeV. This value is consistent with 
the latest preliminary combined results of CDF and DO [18] m t = 171.4 ± 2.1 GeV (y t - 0.94- 0.97). 
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couplings we started with. Therefore, we tune top Yukawa coupling at the splitting scale 
m s and check that obtained values of gi,g2,93 and y t at the electroweak scale are within 
experimental error bars. Also we check that all couplings remain in the perturbative 
regime up to the GUT scale. 

In our restricted parameter space we use RGE to obtain the Higgs boson mass m^. 
The dependence of m h on the splitting scale m s is plotted in Fig. 2. To obtain the plots 

180 
170 
160 
150 

m h , GeV 140 
130 
120 
110 
100 

10 2 10 4 10 6 10 8 10 10 10 12 10 14 
m s , GeV 

Figure 2: The dependence of the Higgs boson mass on the splitting scale m s . Solid lines 
are for cos 2(3 = case, while dotted ones are for cos 2/5 = —1. Top Yukawa coupling 
equals to 0.97 for thin lines and 0.94 for thick lines. 

shown in Fig. 2, we use initial conditions (43) for A and k, and vary the value of (3. In 
this analysis we take into account the experimental uncertainties in the determination of 
the top quark mass. 

The values of the tree level Higgs boson mass are generally within the same range 
as in the minimal split supersymmetry [11, 19]. The upper bound on is increased in 
comparison with MSSM. We have also found that the dependence of the Higgs boson 
mass rrih on k is within 1% in the whole perturbative range of k, while rrih increases 
from 144.6 GeV at A = 0.0 to 160.0 GeV at A = 0.7 (other parameters in this case are 
k = -0.5, cos 2(3 = 0, m s = 10 9 GeV). 

Below we use the two reference points for the relevant parameters of the model, which 
are presented in Table 1. Both of them correspond to rrih — 149 GeV. 





v s 


v P 


h 


A 2 




V- 


(1) 


166 


310 


83.6 


-40.0 


-1.1 


69.6 


(2) 


110 


-310 


-78.3 


-24.8 


-1.0 


-65.2 



Table 1: Two examples of parameters. All values are in GeV. 



y t = 0.97 cos2p = 
y t = 0.97 cos 2p = -l 

y t = 0.94 cos 2p = 
y t = 0.94 cos 2p =-1 
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3 Electroweak phase transition 



For successful electroweak baryogenesis the phase transition must be strongly first order, 
so that the condition v c /T c > 1.1 should be satisfied. Now we turn to the analysis of 
EWPT in our model within the parameter space outlined above. The standard analytical 
way to deal with this problem is based on the methods of finite temperature field theory 
(see, e.g., Ref. [20] and references therein). 

Finite temperature one-loop effective potential for the Higgs and singlet scalar fields 
reads as follows, 

V T (v,vs,v P ) = V tree (v,v s ,v P ) +V {1 \v,v s ,vp) + V^\v,v s ,vp). (44) 

Here the first and the second terms are the tree level part of the potential (31) and 1- 
loop contribution (39), respectively. The third term is the one-loop contribution at finite 
temperature [22], 

= X) /i fa, T) , ( 45 ) 

i 

with 

T 4 f°° / \ 

f t (m i: T) = (±) — J dxx 2 \og (lTe-V^KW), (46 ) 

where we use the same notations as in eq. (39). We found that in numerical calculations 
it is convenient to use the following approximation 

7r 2 T 4 m 2 T 2 m 3 T m 4 m 2 m 
/B(m ' T) = -^ + ^-l2V-64^ ln ^' fOT T <2 - 2 

c . „ s 7tt 2 T 4 m 2 T 2 m 4 , m 2 r m 
^ (m ' T) = -^2b- + ^ + 64^ ln ^' for T <L6 ' 
where In cp ~ 5.41 and In cp ~ 2.64. The approximation 

, / ^ ( rn \ 3 /2 , m ( 15T\ 

is valid for bosons at m/T > 3.0 and for fermions at m/T > 2.2. We adopt linear inter- 
polation between high- and low-temperature regimes. We define the critical temperature 
T c as a temperature at which the first bubbles of true vacuum begin to nucleate. It takes 
place when S3 (T c ) /T c ~ 130 — 140 [16] where S3 (T) is the free energy of critical bubble. 
The critical bubble is a saddle point of the free energy functional, which for spherically 
symmetric configurations reads 



POO 

S 3 (T) = 4tt / rfrr 2 
Jo 



fdh\ 2 fds\ 2 fdp\ 2 T/ o 
+ h- + h- +v T (h,s,p) 



\dr J \dr J \dr J 



Here the effective potential V T is defined by eq. (44), and shifted in such a way as to obey 
Vt (0, S s , P s ) = 0, where S s and P s are the values of the scalar and pseudoscalar fields in 
the symmetric phase. The bubble profile satisfies the equations of motion 

. , d 2 h 2dh dV T . . d 2 S 2dS dV T 



11 



E P (r) = 

and obeys the boundary conditions 



d 2 P 
dr 2 



2dP _dV T 
r dr dP 



(ft(r),S(r),P(r)) 



\dr dr dr J 



=0 = (0,0,0). 



(47) 



To find the critical bubble profile numerically we use the method, originally proposed in 
Ref. [23] for MSSM and further developed in Ref. [24]. Namely, we seek for the absolute 
minimum of the following functional, 

poo 

F [h(r), S(r), P(r)} = An drr 2 [E 2 h (r) + E 2 S (r) + E 2 P (r)] , 

Jo 

over the space of functions satisfying the boundary conditions (47). 

We present the results of numerical analysis for parameters given in Table 1. In fact, 
these values were chosen to have the first order phase transition. In Table 2 we show the 
following quantities: the critical temperature T c , the critical values of the Higgs, scalar 
and pseudoscalar fields, v c , S c and P c , respectively, and also the values of the scalar and 
pseudoscalar fields, S s and P s , in the symmetric phase at the critical temperature. The 





T 

J - c 


Vc 




Pc 


S s 


Ps 


Ss/T 


St (v c ) 


Mo) 


(1) 


88 


218.4 


166.4 


308.4 


350 


68 


130 


0.13 


0.67 


(2) 


96 


215 


110.8 


-307.6 


326.5 


-70.4 


133 


0.18 


0.59 



Table 2: Examples of parameters of the strongly first order phase transition. All dimen- 
sionful values are in GeV. 



critical bubble profile (i.e. the dependence of the scalar fields on radial coordinate) for 
the set (1) of the parameters in Table 1 is presented in Figure 3. The points in Table 2 
correspond to the strongly (v c /T c > 1) first order phase transition. 

Let us comment briefly on the applicability of the perturbation theory for the calcu- 
lation of the effective potential Vr near the critical temperature. Generally, the pertur- 
bation theory at finite temperature is plagued by IR divergences coming from the gauge 
and scalar sectors 6 . The IR divergences in the gauge sector invalidate the perturbative 
calculations at v — 0. However, the small magnetic mass (of order of g 2 T '/ (47r)) of the 
gauge bosons is generated nonperturbatively, and it serves the infrared cutoff for the per- 
turbation theory. It was shown in Ref. [26] that the effect of the magnetic mass renders 
the one-loop effective potential a good approximation even near v — 0. This effect comes 
from the interaction of gauge bosons and in our model it is the same as in SM. As for 
scalar sector, the dimensionless expansion parameter St for unresummed effective poten- 
tial can be estimated as St(v) = ST 2 /6m T (v) [27], where mr(v) is the smallest thermal 
mass of scalar particle (we do not show explicitly the dependence of St on vev's of sin- 
glet scalar fields) and S is the corresponding coupling constant, which, for our parameter 

6 Fermions have no Matsubara zero modes, so their contributions are finite in the infrared region. 
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Figure 3: The critical bubble profile for the set (1) in Table 2. 

choice, is equal to A ~ 0.36 (cf. eq. (35)). We have calculated the thermal mass and 
the expansion parameter at the critical temperature T c for both minima of the effective 
potential Vt and present the results in the last two columns of Table 2. One concludes 
that the perturbation theory is applicable for the points in Table 2. 




4 The baryon asymmetry 

In this section we estimate the amount of the baryon asymmetry produced during the 
phase transition. There are several approaches to this problem (see, e.g., [28], [29], [30] and 
references therein). We have chosen to use the diffusion approximation for the description 
of the particle densities in the presence of moving bubble wall and WKB approximation 
for the CP-violating sources. We take into account CP-violation coming from chargino 
sector only. For the sake of simplicity, we ignore the singlino sector contribution to the 
diffusion equations leaving this for future study. 

According to the semiclassical picture, particles and antiparticles have different disper- 
sion relations in CP- violating background of the bubble wall [31]. Initially, the asymmetry 
in their densities emerges in the chargino sector and then, due to interactions and dif- 
fusion processes, it is transmitted into the densities of other particle species including 
left-handed fermions and, finally, to the baryon asymmetry. 

We derive the diffusion equations along the lines of Ref. [32]. We describe the behaviour 
of the z-th type of particles by the currents = fn^jij, where is the particle density 

and ji is the current density. In the diffusion approximation the current densities have 
the form j% — —DiVrii, where Di are diffusion constants, so diffusion equations can be 
written in the form 

hi - AV 2 ^ = Ii . (48) 
Here the dot denotes the time derivative, 7j are source terms to be described below. To 



13 



the linear order in chemical potentials /ij, one has 

f d 3 k r „ /rn . , /rn kAmi/T)T 2 lA . 

^ = <7* y [/(£*, ^) - f(E h -m)] = %K l/ 6 > fu, (49) 

where ^ is the number of degrees of freedom, f(E i ,iJL i ) is boson (or fermion) distribution 
function, factor k^rrii/T) is equal to 2 (1) for one boson (fermion) massless degree of 
freedom (m 8 = 0) and is exponentially suppressed in the limit rrii/T — > oo. 

First, we neglect CP- violating reactions and weak sphalerons. The CP-conserving 
part of the source terms can be written (again, to linear order in chemical potentials) as 
follows, 

r k 

where the first sum runs over all reactions with the rate T[, in which the i-th type of par- 
ticles participates, and the second sum runs over all types of particle species participating 
in the r-th reaction. As an example, for the process A + 2B — > C + 3D with rate T the 
corresponding term for particle A is I a 3 —Ta {^a + 2/i# — /ic — 3//d)- 

We neglect all Yukawa interactions except for the top quark. Then, leptons do not 
participate in the processes of diffusion of the particle-antiparticle asymmetry. Moreover, 
strong sphalerons, induced by the operator Ytl=i ^Ri u LidmdLi, where index % runs over 
three generations of quarks, are the only processes which produce light quarks; therefore, 
one algebraically constrains the quark densities as follows, 

n Qi = n Q2 = - 2n u R = -2n dR = -2n SR = -2n CR = -2n B , (51) 

where n,Q i = n ULi + n du , (i = 1,2,3) are the densities of left doublets and ng = n^. 
The conservation of the baryon number, along with aforementioned constraints, results 
in tib + tit + riQ = 0, where we use the standard notations tit = n tR and uq = uq 3 . 

Since we neglect the singlet sector contribution, we only write down 7 the source terms 
for the Higgs doublet H' = (H + , H°) T (we use prime here to save the notation H for later 

use), higgsinos H u = (^H+,H^j , H d = (H° d ,H d ^ and heavy quarks Q = (t L ,b L ) T , 
tR and bft . The gauge interactions are assumed to be in equilibrium, so the chemical 
potentials of the gauge bosons are equal to zero. Therefore, the chemical potentials of the 
particles in SU(2)l doublets are equal and, for our purposes, it is sufficient to consider 
only the sums of densities uq = n tL + rib L , uh' = n H + + n H o, rifj u = + n^ and 
n H — nfjo+nfj-. We also make usual assumptions, that the chemical potentials of winos 

« d d 

and bino are zero, being driven by the high rate of helicity flipping interactions induced 
by their Majorana mass terms. 

The following interactions are considered in the diffusion equations: 

(i) strong sphalerons with the rate T ss = 6/t ss |o;fT, where k ss = 0(1) [32], 

(ii) top quark Yukawa interactions y t QtH with the rate Ty, 

(iii) interactions due to the term fiH u eH d with the rate T^, 

(iv) top quark mass effects with the rate T m , 

7 We work in basis of the fields in the unbroken phase. 
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(v) Higgs boson self interactions with the rate Tjj, 

(vi) Higgs-higgsino-gaugino interactions with the rates T u and T d for up and down 
higgsino doublets, respectively. 

We neglect the curvature of the bubble wall and suppose that the main part of the 
baryon asymmetry is produced when the bubble expands already in the stationary regime. 
It means that near the wall all space-time dependent quantities should depend only on 
the combination z + v w t, where v w is the bubble wall velocity and z is a coordinate along 
the axis, perpendicular to the wall, whose positive direction points to the broken phase. 
In particular, n» (f, t) — rij (z + v w t) and hence hi = v w rii. 

With above assumptions one writes the diffusion equations as follows 

v w n' Q -D q riQ = -(fiQ - fi T - [i H ,)T Y - (fi Q - Mt) T m (52) 

-6 (2/ig - fx T + 9/i B ) T ss , 
v w n' T - D q riT = {hq - [It - Hw) Ty + (/jLq - /jLt) T m (53) 
+3 (2/iQ - ht + 9/ii?) T ss , 
v w n' H , - D h n" H , = - (fj, H , -Hq + fJr) T Y - {hh< - HfiJ T u (54) 

- {ji H , + fj,jj d ) T d - Hh'Th, 
v w n'jj u - D h n" 6u = - (nfj u - fi H ,) T u - {ji ku + //^J T M + S u , (55) 
v w n' Ad - D h n" Ad = - {nH d + m>) v d - {vh u + VhJ T M + S d , (56) 

where it is also assumed that the diffusion constants D q are the same for all quarks 
and D h for all higgses and higgsinos. Here we include CP-violating sources S u and S d , 
which we describe later. Neglecting moderate RG effects one estimates the ratio T u /T d = 
tan 2 j3. We consider the limit of large tan /3, when T u ^> T d , and assume the interactions 
corresponding to the rate T u to be in equilibrium, which implies the constraint hh'—^h = 
0. Introducing the notation = T^ + r^ and excluding the terms with T u from eqs. (54), 
(55), one can write down the following linear combinations 

I H'+ I H U = ~ (VH> ~ HQ + fir) Ty - (hh> + HH d ) I> - /J-H'^H + S d , (57) 

Iu d = -^H d + HH')r fi + S u . (58) 
Upon defining the densities 

n h = n H > + njj + rift n H = n w + n k - 



15 



we obtain the following set of diffusion equations 



V w n Q — DgTl'q ~ 


-r y 


n Q _ 
k Q 


n T n H + n h 
k T k Hl 


_ r 

1 m 


>Q 


kx 


-6r ss 


. k Q 


riT 
kx 


, o n Q + n T - 
k B 










v w n' T = D q n'^ + Ty 


~UQ _ 

_k Q 


n T n H + n h 
k T k Hl 


+ r m 




riT 
kx 



(59) 



(60) 



v w n h 



+3r 

D h n" H + IV 



AX + r y 

n H + n h 



k 



kQ k T 

r nQ _nr _ 
.kQ k T 

n Q n T 
kQ k T 

rv + s u + Sd, 



k B 
n H + n h 

k Hl 
n H + n h ~ 



r n H + n h 

J- H — ; r o u — b d , 



k 



k Hl 
n H | n h 
k H k h 



(61) 
(62) 



where 



2(k H f + kfr )kfr 

k Hl = 2(k HI + k &u ), k H = ± h " ' 



k H ' 



fx TT 



k h 



2(k w + kjj n )kfj d 
k ii d + k H> + k &u ' 



We note in passing that in the case tan/5 ~ 1, the assumption, that Y u and Yd are 
in equilibrium leads to a single diffusion equation for n#, which is sourced only by the 
difference S u — Sd- As we will discuss later, in the semiclassical limit S u = Sd, hence 
no baryon asymmetry is produced, however, other approaches to the calculation of the 
CP-violating sources, see, e.g., [29,34,35]) can lead to non-zero difference S u — Sd- 

Making standard assumption that strong sphaleron and Yukawa interactions are in 
equilibrium and follow procedure described in Ref. [29], we find the relations 



k T (2k B + 9k, 



Q) 



(n H + n h ) , n Q 



k Hl (9k Q + 9k T + k B ) 
and the equations for densities n B and 



k Q (9k T - fcg) 
k Hl (9k Q + 9k T + k B , 



A 



= D 



n 



n 



n H 
n h 



+ S, 



(n H + n h ) (63) 



(64) 



where 



Q = 
where f ,■ 



A 

r m + r ft + 2Tfi 
G -r D - F 



1 

F 
F+G 



F+G 
1 



D 



A + A 



A 



D q 
D q + D h 



T m + T h + 2^r, 



D q , D h 



' ' D h and 



o^p (S u + So) 
g+f — Sd) 



(65) 



G+F 1 " ^9 f+G^1-> ^ n F+G J 

F = 9k Q k T + k Q k B + 4k T k B , G = k Hl (9k Q + 9k T + k B ) 
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Using (63), one finds the density of left-handed fermions 



n L eft = n Ql + n Q2 + n Q3 = 4n T + 5n Q 

_ 5k Q k B + 8k T k B - 9k Q k T 
~ k Hl (k B + 9k Q + 9k T ) 



(n h + n H ) = A- (n h + n H )- (66) 



For our set of statistical factors (k,Q = 6, kr — 3, k B = 3, kw = 4, kjj u = k H - d = 2 in the 
massless case), which is the same as in SM (except for the higgsino sector), the value of 
constant A equals zero. This is well known Standard Model suppression [36]. As shown in 
Ref. [14] for similar situation, the corrections due to subleading 0(l/r ss ) effects are small, 
while radiative corrections to the thermal masses of particles give a sizable contribution 
mainly due to the top Yukawa coupling, 

The space-time distribution of left-handed fermions acts as a source for baryon asymmetry 
[29] , which obeys the following equation 

v w n' B {z) = -0(-z) [n F T ws n Left {z) + Kn B ] , (68) 

where Y ws is the weak sphaleron rate (T ws = 6/c ws a^T, where k ws = 20 ± 2 [37]), 1Z is 
the relaxation coefficient [28,38] (TZ = ^npT ws ), and rip = 3 is the number of families. 
Here theta-function implies that the weak sphalerons are active only in the unbroken 
phase. The solution of this equation can be obtained analytically, and the resulting 
baryon asymmetry takes the form 

nB = _!^£ f° dzn Left (z)e^ . (69) 

In the WKB approximation, the source terms 8 S u and Sd are derived from chargino 
dispersion relations. The mass matrix of charginos reads as follows, 

M ch =[ 2 ), 70 

where we introduce the notation p,(z) = ii + k(S(z) +iP(z)) The expressions for 

the sources (up to overall factor) come from Ref. [8], where it has been shown that, in 
fact, the CP- violating sources for higgsinos H u and are equal, so only the equation for 
density rih is sourced. The expressions for the sources read 



S * = ^T^{^ A ) \- 6 ' - V sin2 a + 6 ' sin2 b 1 '+ 



(B) (ml)' [9' - i sin 2 a + 5' sin 2 b] }', (71) 



8 Here we include CP-violating sources only for charginos. We make the standard assumption that 
the contribution of neutralinos produce a correction by a factor of order unity to the final result for the 
baryon asymmetry (in MSSM this factor is about 1.5). 
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where 

(A\ = /\P 

2E 2 I + \ 2E A 

and E 2 = p 2 + p z is the component of the spatial momentum perpendicular to the 
wall, and m~ h is the mass eigenstate of the matrix (70), which becomes pure higgsino in 
the unbroken phase. The thermal averages ( ) + and ( ) are performed with fermion 
distribution function 

um Jd? P (df/dE)A(P) df _ eg 

{Am ~ f d *p(df/dE)\ m=0 ' dE~ ^(e^ + 1) 2 {76) 

for massive (+) and massless (0) cases, respectively. The other quantities in eq. (71) come 
from the diagonalization of chargino mass matrix [8], 

sin 2 a = 2|A| 2 /A(A + A), A = h (g u M 2 + g d Jj) /V2, A = (A 2 + A\A\ 2 ) 1/2 , (74) 
7 = argA 7 / sin 2 a = 2Im(AM')/A(A + A), A = M 2 - |/i| 2 + {~g 2 d - ~g 2 u ) h 2 (75) 

The expressions for sin 2 b, 5 and 5' sin 2 b are obtained from the ones for sin 2 a, 7 and 
7' cos 2 a by replacement a — > b, g u — > g d and 7 — > —5. 

To calculate CP-violating sources, we need the bubble wall profile. The standard 
way to find it is to solve the equations of motion, derived from the effective potential 
Vr{h, S, P) at the critical temperature T c (for recent calculations of the bubble wall profiles 
in (N)MSSM see, e.g., Refs. [24]). For our purposes we approximate its form by the 
standard kink-like solution 



h(z) = ^v c 

S(z) \ = ( S c \ _ If Av s 
P(z) J \PJ 2\Av P 







^1 — tanh 









1 + tanh 



2z 

all- — 

L>,„ 



(76) 
(77) 



where S c , P c are the critical values of the singlet scalar fields, Av s = S c — S s and Av P = 
P c — P s (see also Table 2), L w is the bubble wall width, which is taken to be the same 
for all scalar fields, and the constant a is taken to be 1.5 as in Ref. [29]. Here we have 
assumed that at z > L w j (2a) gauge symmetry is broken (the inner region of the bubble) 
and at z < symmetry is restored (outer region). 

For numerical calculations we take the diffusion constants to be D q ~ 6/T for quarks 
and Dh ~ 110/T for Higgs sector [33]; the width and the velocity of the bubble wall are 
taken as L w = 8.0/T and v w = 0.1, respectively. The rates 

T H = 0.0036T 9 (z - 0.5L W ) , T m = 0.05T 9 (z - 0.5L W ) (78) 

are the same as in Ref. [8] and = OAT [7]. 

Below we present the results for the baryon-to-entropy ratio A = ub/s with entropy 
density s given by s = 2ir 2 g e ffT 3 /4:5, where g e ff is the effective number of relativistic 
degrees of freedom at the critical temperature. We plot the ratio A/A , where A = 
8.7 • 10~ n , which corresponds to n B /n 7 = 6.5 ■ 10~ 10 (cf. eq. (1)). In Fig. 4 we show 
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M 2 , GeV 

Figure 4: Plot of A/A as a function of M 2 for the set (2) of parameters presented in 
Table 2. 

the dependence of the baryon asymmetry on the gaugino mass M 2 for the set (2) of the 
parameters in the Tables 1 and 2. We found that in contrast to MSSM, the baryon 
asymmetry does not vanish in the limit of infinite M 2 , where one of the chargino species 
decouples. This fact may be attributed to the strong z-dependence of the effective jl 
parameter. We have checked numerically that in the limits M 2 — > +oo and M 2 — > — oo 
the results for the baryon asymmetry coincide as it should be, since in both cases charged 
gaugino (wino) is very heavy and does not affect the low energy physics. Similarly, 
at \M 2 \ << the amount of baryon asymmetry is not very sensitive to M 2 and is 
determined by jl(z) and off-diagonal entries of the mass matrix (70). In both regions 
(M 2 — > oo and M 2 — ► 0) wino admixture in the relevant higgsino-like mass eigenstate 
is small, so our approximation where we neglect mixing in diffusion equations can be 
justified. The region of large higgsino-wino mixing (IM2I ~ ~ 200 GeV in Figure 4) 
should be considered more accurately if one wants to know the exact amount of baryon 
asymmetry produced in this case (see [30] and references therein). 

To summarize this section, we demonstrated that the model is indeed capable of 
producing sufficient amount of the baryon asymmetry. There are certainly values of 
parameters for which the electroweak phase phase transition is sufficiently strongly first 
order, and electroweak baryogenesis is sufficiently efficient. In this analysis we restricted 
ourselves to rather small range of parameters. It remains to be understood how generic 
is successive electroweak baryogenesis in the whole parameter space. 

5 Electric dipole moments 

One of the consequences of the presence of additional sources of CP-violation is their 
contribution to the electron and neutron EDMs. Like in the minimal split SUSY model, 
there are three types of contributions to EDM of a fermion (lepton or quark), related 
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to the exchange of h -y, W + W~ and h Z bosons (see Refs. [39], [40], [41], [42]). The 
corresponding two-loop Feynman diagrams for SM fermion / (charged lepton or quark) 
are given in Fig. 5. It was shown in Ref. [42] that the third contribution is suppressed 




Figure 5: Feynman diagrams contributing to the fermion EDM in split SUSY. 



in the case of electron by a factor (1 — 4 sin 2 &w), where 9w is the weak mixing angle, 
while it is important for the neutron EDM. In the minimal split SUSY model there is 
only one CP-violating phase, associated with the invariant 0i = arg(g u * g d * M 2 fl) where 
jl = fx+K (vs + ivp) I In the class of models presented here, among parameters directly 
governing the couplings contributing to the diagrams in Fig. 5, two additional independent 
phases appear, which correspond to the phase invariants 02 = arg i^k* \ u \d(fi*)~ 2 ) and 
03 = arg(A u A^ d ). 

After diagonalization of chargino and neutralino mass matrices, the relevant part of 
the interaction Lagrangian becomes 

9 



-Ant = -^Xtl" {GfjPR + G%P L ) xp, 

9 - 



(79) 



+ gxt-f {C*Pr + C^P L ) x jW+ + ^=xt (D*P R + D%P L ) x +h + h.c. 

where the matrices are 

G% = Vacw+V^ + V l2 c hi V 2 ), -G% = U lX c w -U\ 3 + U l2 c h -U\ 3 , 
gD% = ~g* a V t2 U n + ~g* d V tl U j2 , gD L = (gD R ) ] , 
with i,j = 1,2 and c w ± = ± cos 2 9 W , c h + h - = ±(1/2 — sin 2 9\y), 



r L - 



-V a N. 



P- 



~ u * lNj2 " Tf^' 3 ' 



(80) 
(81) 

(82) 



with i — 1,2 and j = 1, . . . , 5. Here the unitary matrices V, U and N diagonalize the 
mass matrix of charginos (70) and the mass matrix of neutralinos 



M n = 
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Mi 





~\9> 









M 2 


\~9dV 








\9dV 





-jl 






\g u v 















T2 XdV 








7S A ^ 



\ 



(83) 



j=X u v -V2k(v s + iv P ) J 
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m „+) and N*MN* 

A2 



diag (m x o , m x o , m x o , m x o , m x o ) . The 



so that {/"".Mc^ = diag(m 
EDM of electron or light quark has the form 

d f 

where expressions for the corresponding contributions can be found in Ref. [42]. We 
modified them appropriately in our case, because of different number of neutralinos. The 
neutron dipole moment can be expressed in terms of quark EDMs by the relation [43, 44] 



4 = (1±0.5) 



4 , 
-d d 



—d, 
3 



(84) 



(m„ + m d )(225MeV) 3 _ 
where /„- = 92 MeV and m u /m d = 0.55. 

For numerical calculations we use two sets of parameters (see Table 1) and randomly 
scan over the following parameter space, < Mi,M2 < 1000 GeV. Also we take the 
coupling k to be complex, k = \k\e ir ^ k , (0 < <pk < tt) to include the contribution of the 
phase invariant <p2- The results for the electron and neutron EDMs as functions of the 
mass of the lightest chargino are presented in Figs. 6 and 7, respectively. We have taken 
into account the experimental bound on the mass of the lightest chargino m x + > 104 GeV 
[45]. 



le-27 



ld e /el, cm 



le-28 



1.6 10" 27 cm 





1.6 10" 27 cm 


le-27 






ld e /el, cm 












le-28 







100 120 140 160 180 200 220 100 120 140 160 180 200 220 

VRy+, GeV m^+, GeV 

Figure 6: The EDM of electron as a function of the mass of the lightest chargino m x + for 
set 1 (left) and set 2 (right) of parameters in Table 1. Horizontal dotted line represents 
the experimental bound \d e \ < 1.6 • 10 _27 e cm. Green (gray) points represent the values 
of EDM for real k. 

Horizontal lines show the present experimental limit on the EDM of electron \d e \ < 
1.6 • 10" 27 e cm at 90% CL [46] and neutron \d n \ < 3.0 ■ 10~ 26 e cm at 90% CL [47]. One 
observes that generally for the cosmologically favorable models, the predictions for EDMs 
are within one-two orders of magnitude below the present experimental bound. Hence, 
our solution of the baryon asymmetry problem can be indirectly probed by the future 
experiments aimed at EDM searches. 

6 Dark matter candidates 



In this section we explore dark matter in the model. In the minimal split SUSY this 
issue has been already investigated in Refs. [12], [39], [48]. In the first place, let us note 
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Figure 7: The EDM of neutron as a function of the mass of the lightest chargino m x + 
for the same sets of parameters as in Fig. 6. Horizontal line represents the experimental 
bound \d n \ < 3.0 ■ 10~ 26 e cm. Green (grey) points represent the values of EDM for real k. 



that a generalization of R-parity can be introduced in the nonminimal split SUSY: with 
respect to this R-parity all new fermionic fields are odd, while new (pseudo)scalars are 
even. Hence the lightest new fermion is the lightest superpartner and it is stable. 

The viable dark matter particles should be electrically neutral, so the best candidates 
in our model are neutralinos. We define the LSP neutralino state as 

X = N 51 B + N 52 W + N 53 H U + N 54 H d + iV 55 n. (85) 

We estimate the neutralino relic abundance by using the standard methods [49]. First we 
calculate the freeze-out temperature Tp, 



( m / 45 \ 

^ =1o h^v^^ Mm01 J' (86) 

i ii 

where Xp = m x /Tp, g* is the effective number of degrees of freedom at freeze-out gj ~ 9; 
the thermally-averaged product of neutralino annihilation cross section a and relative 
velocity v of neutralinos (the M0ller cross section (<tv)m0i) is 

1 f°° ( fs\ 

MM01 = 8<T^K/T) L x dMsKS " A < )V7sK ' {-T ) ■ (87) 

Here K\, K 2 are the modified Bessel functions, s is the Mandelstam variable. The relic 
abundance of the lightest neutralino is then given by 

= (1.07xlO»GeV-) / r dx M X ) 1/2 
x Me, \J IF 

To calculate the neutralino abundance, we modify the formulas for the annihilation cross 
sections presented in Ref. [50]. The changes concern the annihilation of neutralinos into 
the Higgs bosons and are due to new Feynman diagrams with singlet scalar and pseu- 
doscalar exchanges in s-channel. The corresponding modifications are presented in Ap- 
pendix B. The total number of neutralinos has been also changed according to our model. 
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Adopting constraints discussed in Sec. 2, we first scan uniformly over the following 
parameter space: \Mi\, |M 2 | < 1000 GeV, \vp\ < 2000 GeV, with vs being in the region 
of correct electroweak vacuum (i.e. the electroweak breaking minimum is the global 
minimum of the potential) and squared mass matrix of the scalar fields being diagonal. 
The numerical results are presented in the left plot in Fig. 8, where we show the region 
in (m x +, m x ) plane favored by WMAP data: each point corresponds to a model in which 
neutralino abundance is within the range 0.094 < QoMh 2 < 0.129 [1]. To check that the 
baryon asymmetry and dark matter problems can be solved simultaneously, we also scan 
uniformly over the region in the parameter space preferred by electroweak baryogenesis 
(cf. Sec. 4); namely, we use \M\\, \M%\ < 1000 GeV, with other parameters corresponding 
to the set (2) in Tables 1 and 2. Points in (Mi, M2) plane, which correspond to correct 
neutralino abundance, are shown on the right plot in Fig. 8. 




Figure 8: Points in (m x +, m x ) plane, i.e. the lightest chargino and the lightest neutralino 
(left) and in (M 2 , Mi) plane (right), which correspond to models with the correct relic 
abundance of singlino (|AT 55 | > 0.5, green (light grey) oblique crosses) and non-singlino 
(red / dark grey crosses) dark matter. 

On both plots green (light grey) crosses correspond to the dark matter particles which 
have considerable admixture of singlino (| JV55 j > 0.5), while the red (dark grey) crosses 
correspond to the mostly bino LSP. The annihilation of DM particles (bino as well as 
singlino) with masses m x ~ 0.5Mh ~ 75 GeV or m x ~ 0.5M^ proceeds resonantly 
via Higgs or Z°-boson exchange, respectively. On the right plot, this light neutralino 
corresponds to the red (dark grey) horizontal lines with Mi < 80 GeV. The most part of 
the parameter space with singlino-like dark matter give relatively light LSP with mass in 
the range 50 — 200 GeV although heavier candidates are not entirely excluded. We have 
found that in this case, considerable admixture of higgsinos is always present (numerically, 
we obtain ( | 1 2 + I./V54I 2 ) ^ ~ 0.4 — 0.8). Singlino dark matter with the mass m x > 
80 GeV annihilates predominantly into W + W~ gauge bosons, while for m x ~ 0.5Mh or 
0.5Mz the main channel is the resonant one, XX ~^ h*(Z *) — > //. One concludes that 
dark matter problem can also be solved in the framework of considered models. 
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7 Discussion and conclusions 



In this paper we proposed a generalization of the minimal split supersymmetry model by 
taking into account the necessity to explain the baryon asymmetry of the Universe. The 
main idea was to start the whole construction of split SUSY not with MSSM but with 
NMSSM. We needed to fine-tune some parameters of the model to obtain the same low 
energy spectrum of particles as in the minimal split SUSY plus additional singlet particles, 
which give rise to all features one needs for successful baryogenesis. We restricted ourselves 
to rather small part of parameter space, in which the lightest Higgs boson does not mix 
with other scalar particles. Like in the minimal version of split supersymmetry model, the 
upper bound on the mass of the lightest Higgs boson shifts to larger values in comparison, 
e.g., with MSSM. 

We have considered the electroweak phase transition and, by exploring one-loop ef- 
fective potential at finite temperature, we found that there is a region in the parameter 
space where the phase transition is strongly first order and the baryon asymmetry is not 
washed out after the phase transition is completed. We have used WKB approximation 
to calculate the value of the baryon asymmetry. We have found that this model is indeed 
capable of producing the right amount of the baryon asymmetry for realistic values of the 
width and velocity of the bubble wall. 

We have investigated the contribution of the CP-violating sources into the electron 
and neutron EDMs. As in the minimal split supersymmetry, their values are close to the 
present experimental limits and can be measured in the future experiments hence giving 
an opportunity to falsify the suggested solution of the baryon asymmetry problem. 

We have explored the dark matter in nonminimal split SUSY. In addition to the bino 
and higgsino dark matter candidates there is a region in the parameter space in which 
singlino can considerably contribute to the LSP state and hence plays the role of dark 
matter. We have found that in the latter case LSP is quite light (m x < 150 GeV) and 
contains considerable admixture of higgsino as well. 

To summarize, the split NMSSM models are capable of solving both baryon asymmetry 
and dark matter problems and can be probed by the next generation of EDM experiments. 
The collider phenomenology of this model is quite similar to one of minimal split SUSY, 
if singlino-neutralino and higgs-singlet mixing is small. In the opposite case there are 
additional signatures of this model resembling ones in non-split NMSSM. We leave the 
study of LHC prospects in probing this model for the future. 
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8 Appendix A 



This appendix contains the set of renormalization group equations for the coupling con- 
stants of the model described by the Lagrangian (31), (32) at energies below m s . One 
derives them by making use of the general formulae given in Ref. [51]. 
The 2-loop RGE for gauge couplings have the following form 

d _ _sl a! 



(^)j^ = 9% + ^- 2 {^B l ^- J2 d?Tr(h«h°) (89) 

^ ' j=l a=u,d,e 

-€ (gl + <9 - df K + sfi) - 4 (M + a 2 ) - d^) , 



where t — In// and Jx is the renormalization scale. Here the GUT convention g\ = (5/3) g' 2 
is used. Below we list the coefficients of /9-functions relevant within the interval m ew < 
p, < m s . Above the splitting scale, m s , the renormalization group equations become the 
same as in the usual NMSSM (see, e.g., Ref. [52]). One has 

( 104 18 44 \ 
25 5 5 \ 
,90) 

H^H- <^frH- '=(§.H- (91) 
dW -{^»' M^H- ""-(I- 1 - )- « 92) 

It is convenient to introduce the following notations 

T H = Tr [3h«h» + 3h«h d + h*h e ] +\{gl + gj) + \ (g 12 + ~g' d 2 ) + A 2 + A 2 , (93) 

T N = 2 (k 2 + k 2 ) . (94) 

One-loop equations for the quark and lepton Yukawa couplings can be read off from Ref. 
[12] with the substitution T — > T H . One obtains 

(An) 2 ^-h u = h u (-3 J2 ^ + - !W + T « ) > (95) 



i=i 



d , j , j I „ x ^ ri o 3 , „,+ , „, 3 



(4tt) 2 = /i d -3 V cf # 2 - -h u ^h u + + T H , (96) 

at ' ^— ' " 9 / 



=i 



eft 

Here the coefficients c"' d ' e are given by 



(4vr) 2 -/> e = \f -3 X] <fr? + -/^ e + . (97) 



i=i 



1Z 3 8\ WJ_3 8\ c e = ^3 
60'4'3;' ^12'4'3y ' V4'4' 
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RGE for gaugino and singlino couplings (see (38) and (36)) are 

d ^ 5 1 X 

(4vr) 2 Jt ~9u = -3~ 9u °d + ~/u - ^9u9 2 d + + MWu (99) 

i=i 

1 1 

+ 2#« A « + ^9uK 2 + 2g d \ u \ d + g u T H , 

3 

(4vr) 2 j t ~g' u = -3g' u £ C'g 2 + \g« + \g'Jj + \g' u g\ + 3^ (100) 
3 _ ~ 1 

+2^ A « + 2^ + Zg' d \u^d + g'u T H, 

(An) 2 j\ u = + 3AX + 2k\ u + \\ u k 2 + \\ u (3g 2 u + g' 2 ) (101) 

3 



+A d (3M« + g'Ju) - 3A„ ^ C[g 2 + X u Th, 



i=l 



where 

C7=( 

V 20 ' 4 7 V ~ V20' 4 : 
The equations for ^ and A^ are the same as (99), (100) and (101) with the only 
replacement u <-> d. 

RGE for singlino Yukawa couplings k and /c are 

(An) 2 jk = 4A; 3 + 2k (\ 2 U + \ 2 d ) - l -\ u \ d K + kT N: (102) 

(4-) 2 jf = ^ + ^ (A 2 + AS) + \k (g 2 u + sr3) + \« K + (103) 

+4/cA M A d - 3k + |^ + kT at- 

RGE for the Higgs quartic coupling is 

(4tt) 2 |A = 12A 2 + 2k 2 + 8«i + A [-9 (f + g^ +6(g 2 u + ~g 2 d ) + 2 + ~g> 2 ) (104) 

+ 9 , H , <&g" 



+4 (A 2 + A 2 ) + 4Tr (3/i ut /i u + 3h d ^h d + /i et /i e ) j } s >? ( 

-4Tr (3 {h^f + 3 (^^) 2 + (h^f) -5{gt + ft) - 2g 2 u g 2 d - {g' 2 + g' 2 ) 2 

-2 (g u g u + gdg'd) - 4 (A 2 + A^) 2 - 4 (\ u g u + \ d g d f - A (\ u g u + \ d g' d ) 2 . 
Finally, RGE for the remaining scalar couplings have the following form 

(An) 2 ^-ki = 8Ajv/ti + Ak\ + 16k 2 . + 6/^A - At]k 2 - 8k 2 (A 2 + A 2 ,) (105) 

+2k, (T h + T N ) - + 3 2 2 ) , 
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(4tt) 2 = 6A/t 2 + 8/ti/t 2 + 24^2 + 4/t 2 Ajv - 4/c/tA M A d (106) 
+2/€ 2 (T H + T N ) - |k 2 ^ + ^ , 

(4tt) 2 — Ajy = 20A 2 V + 144kik 2 + ^ + 2l 4 + 2l 4 - ^ ~ 2 ^ + 4\ N T N (107) 
at 2 

1 Q 

(4vr) 2 -i = 24A^ + Ak 2 2 + -rf + 4^, (108) 

(4tt) 2 ^-r) = 36rj\ N - Y1k x k 2 + 72^ + Ar]T N . (109) 
at 

9 Appendix B 

This Appendix contains the relevant set of formulas describing the cross sections of neu- 
tralino annihilation into the Higgs bosons. Here we adopt the notations of Ref. [50] for 
couplings and matrices. The relevant part of the Lagrangian is 

C = \fl $ (cf^i + cf^ 2 $ 2 ) X° 3 + \ (C hh ^ + C hh ^ 2 ) h\ (110) 

where we use $i, $2 for mass eigenstates of scalar S and pseudoscalar P fields related to 
CP-eigenstates as follows 

S \ _ ( cos6»$i + sin6»$ 2 
P ) ~ \ -sin6»$i + cos6>$ 2 

with a mixing angle 9 and 

c x°x°^i = _ yj2kN* b N* b e~ iB , Cs^ 2 = -iV2kN* 5 N* 5 e~ ie , (111) 

C hmi = (A 2 cos9 + A 1 sm9) + ]- [( Kl - k 2 ) v s cos 9 - («i + k 2 ) v p sin 9} , (112) 
y2 \ / 2 

— k 2 ) v s sin 9 + (ki + k 2 ) vp cos 9} (113) 

For simplicity we assume, as everywhere in this paper, that the Higgs state h does not 
mix with other (pseudo)scalar particles. We present functions Whh, determining the cross 
sections of nonrelativistic scattering XX ~^ hh (see Ref. [50] for definitions of these and 
other related quantities). There are contributions to the cross sections coming from the 
s-channel Higgs (h), scalar (S) and pseudoscalar (P) exchanges and t- and -u-channel 
neutralino exchange as well as from their interference, 

= + (H4) 



For scalar exchanges in the s-channel one has 

2 

hhrnXXr 

(s-4m 2 ), (115) 



Whh 



(/i,*i,* 2 ) _ 1 



rihhrrixxr 



s — m% + iT r m r 
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The scalar - neutralino interference gives 



w 



(/l,$l,$2-X°) 

hh 



2 E Re 



i=l 



E 

T=fe,*l,*2 



QhhrQXXr 

s 

s — ml + iY r m 7 



(116) 



x |-2m x + [2m x (m 2 h - m^ - m 2 x ) + m x o(s - 4m x )]jF(s, m x , m^,m^,m^ )| 



The term corresponding to t- and w-channel neutralino exchange is the same as in 
Ref. [50]. 
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